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We complement and extend our work on fluctuation relations arising in nonequilibrium systems in 
steady states driven by Levy noise [Phys. Rev. E 76, 020101(R) (2006)]. As a concrete example, we 
consider a particle subjected to a drag force and a Levy white noise with tail index fi G (0, 2], and 
calculate the probability distribution of the work done on the particle by the drag force, as well as the 
probability distribution of the work dissipated by the dragged particle in a nonequilibrium steady 
state. For < ^ < 2, both distributions satisfy what we call an anomalous fluctuation relation, 
characterized by positive and negative fluctuations that asymptotically have the same probability. 
For fi = 2, hy contrast, the work and dissipated work distributions satisfy the known conventional 
and extended fluctuation relations, respectively, which are both characterized by positive fluctuations 
that are exponentially more probable than negative fluctuations. The difference between these 
different fluctuation relations is discussed in the context of large deviation theory. Experiments that 
could probe or reveal anomalous fluctuation relations are also discussed. 
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I. INTRODUCTION 

Many works on nonequilibrium systems have been de- 
voted recently to the study of fluctuations around out- 
of-equilibrium steady states. The basis of these studies 
is typically the following: for a given system driven in a 
nonequilibrium steady state, one considers an observable 
A-r integrated over a time r, such as the work done by a 
force acting on the system during a time r, and proceeds 
to calculate its probability distribution P{Ar). Inter- 
estingly, what has been found for many different systems 
and many observables A-r is that the positive fluctuations 
of Ar are exponentially more probable than its negative 
fluctuations, in the sense that 
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where c is some constant that does not depend on r, and 
o(t) is a sublinear correction term in r, which may be 
zero if the probability ratio is exactly exponential in r. 
Examples of observables for which Eq. ([T]) is observed 
include the entropy production of chaotic deterministic 
systems P, [H, Q and stochastic Markov systems [1, d, Q , 
as well as other entropy- and work-like quantities defined 
in the context of particles moving in fluids H, B, el ectri- 
cal circuits [1, [lol. granular media [ll|, [13111^ , and 
turbulent fluids [1^, ^M, , among other systems. 

In early studies of these systems, it was thought that 
the result expressed by Eq. H]) might be a general law 
of nonequilibrium fluctuations, but it is now known that 
this is not the case. Some observables, such as the heat 
absorbed by driven Brownian particle s \vk fisl . [Tot and 
the current of the zero-range process j20l. l2ll|. satisfy a 
more general relation of the form 
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where /(a) is in general a nonlinear function of a, which 
does not depend on r. The difference between Eqs. ([T]) 
an d dp serves, following the work of van Zon and Cohen 
[1, Ell, ^ ^ basis for classifying nonequilibrium fluctua- 
tions. Observables that comply with Eq. ([T]) are said to 
satisfy a conventional fluctuation relation, whereas those 
satisfying the more general Eq. ^ are said to satisfy an 
extended fluctuation relation. 

Our goal in this paper is to revisit this classification of 
nonequilibrium fluctuations. Based on the fact that the 
existence of a conventional or extended fluctuation rela- 
tion is essentially equivalent to the existence of a prob- 
ability distribution having a large deviation form [2^ . 
we have constructed in [23[ a model of a driven non- 
equilibrium system for which the mechanical work WV 
done over a time r by the driving force satisfies nei- 
ther the conventional nor the extended fluctuation re- 
lation because the probability distribution P{Wr) of Wr 
fails to have the form of a large deviation probability. 
The probability distribution P(WV) is nevertheless well 
defined, and can be used to define the probability ra- 
tio P{Wr) / P{~Wr), which has in this case a power-law 
rather than an exponential form in r. We have called 
this property of the probability ratio an anomalous fluc- 
tuation relation, following a terminology used in studies 
of Levy-type noise, and have proposed some experiments 
with which one could physically "realize" or "test" this 
type of fluctuation relation. 

Here we complete our study of anomalous fluctuation 
relations initiated in [23t by discussing in more detail the 
fluctuations of Wr for the model studied in that paper, 
and by considering the fluctuations of an additional ob- 
servable, called dissipated work. These two points are 
the subject of Sees. IIIII and llVi respectively. We also dis- 
cuss in Sec. |V] the difference between conventional and 
extended fluctuation relations, on the one hand, and the 
anomalous fluctuation relation, on the other hand, from 
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the general point of view of large deviation theory. The 
main result discussed in that section follows the observa- 
tion that fluctuation relations of the conventional and ex- 
tended types are equivalent to having probabilities of the 
large deviation kind, and that other types of fluctuation 
relations must arise whenever the probability distribution 
of an observable does not have the large deviation form. 
The anomalous fluctuation relation that we discuss here 
is but one example of fluctuation relations, which arise 
by replacing Gaussian white noise as the source of noise 
in Langevin equations by Levy white noise or, more gen- 
erally, by replacing noise sources having finite moments 
by noise sources having infinite moments. 

The relationship between fluctuation relations and 
large deviation theory was noted in the original deriva- 
tions of the fluctuation theorem for the entropy produc- 
tion [l], 0, Q (see also [241 ). and is explicitly discussed 
in our previous paper |23| . as well as in a recent review 
paper written by one of us . Here we continue these 
discussions by studying some conditions under which the 
anomalous fluctuation relation is expected to arise. We 
end the paper in Sec. IVII with various remarks related 
to the physical interpretation of our results, the nature 
of Levy noise, as well as future work aimed at extending 
and experimentally verifying our results. 

II. MODEL 

We consider in this paper a Brownian particle sub- 
jected to three different forces: a linear restoring force 
arising from a particle-confining harmonic potential mov- 
ing at a constant speed, a friction force, and a random 
force or noise. The Langevin equation modelling the ef- 
fect of these forces on the Brownian particle can be found 
in the work of van Zon and Cohen [1], which is itself 
based on the experimental work of Wang et al. ■ Here 
we study the overdamped version of that equation, given 

by 

ax = -K[x{t)~vt\+£,{t). (3) 

In this equation, x{t) denotes the position (in the labo- 
ratory frame) of the Brownian particle at time t, v is the 
velocity with which the harmonic potential moves, a is 
the friction coefhcient, k is the strength of the harmonic 
potential, and ^{t) is the random force. 

In previous studies of the model defined above, as in 
most studies of Langevin equations, the random force ^{t) 
is assumed to be a Gaussian white noise characterized by 
its zero mean {C{t)) = and its autocorrelation function 
{£.{t)£.{t')) = r(5(i - t'), where F is the noise power. We 
depart from this assumption here by taking ^(t) to be 
a symmetric Levy white noise, defined by the following 
characteristic functional: 

G^[k] = J V[^] P[^] e'/'=(*)«W'^* = exp (^-b J |fc(t)|^di^ 

(4) 



with & > and < ^ < 2 [H, [H, [l^. The first integral 
in this expression is the path integral defining the char- 
acteristic function of the noise process ^ (t) , whereas the 
second expression on the right-hand side is the expres- 
sion that defines ^{t) as being a Levy white noise with 
strength b and index /i; that is, an uncorrelated noise 
distributed according to a symmetric Levy distribution 
with scale parameter b and tail index /x.^ The case of 
Gaussian white noise considered in Q is recovered from 
this characteristic function by choosing /x = 2, in which 
case ^{t) = and {^{t)^{t')) = TS{t - t') with L = 2& as 
the noise power. 

It is important to note that the use of Levy noise in 
Langevin equations is often seen as problematic, since the 
mean of ^(i) diverges for ^ G (0, 1), and its noise power 
F diverges for aU ^l € (0, 2) [2l,[2i,l23|. The fact is, how- 
ever, that these divergences do not lead to any physical 
pathologies; they are merely the sign that the concept 
of mean and noise power (viz., variance) are ill-defined 
or inapplicable for Levy noise. This point has been dis- 
cussed in the literature; see, e.g., [1^. In particular, the 
fact that ^(t) has an infinite mean for ^ £ (0, 1) does not 
imply physically that the energy supplied to the Brow- 
nian particle by the random force £,{t) is infinite. The 
increments of ^(t) are necessarily always finite, so that 
the energy supplied by ^(t) is also always finite. 

As for the divergence of the power of ^{t) for /i G (0, 2), 
it is true that it implies that there can be no fiuctuation- 
dissipation relation relating the friction coefficient with 
the noise power. However, if we view the friction and the 
noise as being independent, that is, as arising physically 
from different physical mechanisms, then there is no need 
for a fluctuation-dissipation relation between the friction 
and the noise. This situation is physically possible: one 
can imagine, for example, that the friction force in the 
Langevin equation arises from a solid-solid contact, while 
the noise is imposed externally, say, using a computer- 
generated Levy noise. This situation will be discussed 
more concretely in Sec. IVII 

For now we will leave these considerations aside to fo- 
cus on the nonequilibrium steady state generated by the 
Langevin equation defined in Eq. and to study the 
fluctuations of the work done on the Brownian particle 
in such a state by the moving potential. This is done in 
the next section. The fluctuations of the dissipated work 
are studied in Sec. IIVI 



III. WORK FLUCTUATIONS 

The total work done on the Brownian particle de- 
scribed by Eq. ([3]) is the sum of two contributions: the 
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mechanical work done on the particle by the random force 
£,{t) [13 J which is explicitly considered here as an inde- 
pendent external force, and the mechanical work done on 
the particle by the moving harmonic potential Q. We 
study in this section the latter quantity which is given 
by the integral 



where 



Wr 



KV / [x{t) — vt] dt 



(5) 



for a given time interval [0,t]. Thus Wr is the mechani- 
cal work done on the Brownian particle by the harmonic 
potential over a time r. For convenience, we rewrite this 
quantity as 



Wr = -KV / y{t)dt 



(6) 



in terms of the position y{t) = x{t) — vt of the particle in 
the frame of the moving potential (comoving frame) Q. 
In terms of y{t), the Langevin equation of Eq. jS]) reads 



y{t)^--y{t)-v + -m, 

Tr a 



(7) 



where we have defined = a/K, the characteristic re- 
laxation time of the particle in the potential. 

Our aim in the next subsections is to calculate the 
probability distribution of Wr for various values of ^, 
and to discuss the properties of this distribution in the 
light of fluctuation relations. 



A. General distribution 

The calculation of the distribution P{Wr) is done in 
three steps. First, note that the characteristic function 
of Wr can be obtained from the characteristic functional 
Gy[k] of y{t), defined as 



Gy[k] = ^cxp ^i J k{s) y{s) ds 
by choosing the test function 



-QKV, 



0, 



< s < r 

S > T. 



From the definition of the work, given in Eq. 
indeed have 



(8) 



(9) 



we 



Gv^.(g) = (e"'^'-) 



exp 



-iqnv I y{t) dt 



--Gy[k] 

(10) 

In the second step, we express the characteristic func- 
tion of y{t) in terms of the characteristic function of the 
noise process ^ [t) . The Langevin equation for y[t) , shown 
in ([7]), is fincar in y{t), so we can use for this purpose a 
general result of Caceres and Budini [2^, [2^ , given by 

Gy[k] ^ exp (^iroyo - iv r{l) dl^ G^[r/a], (11) 



r{l) 



(12) 



ro = r{0), yo = y{0) = x{0), and G^ is the characteristic 
function of the noise f(<). Substituting the expression of 
G^, given in Eq. we then obtain 



Gwr (q) = cxp iroyo 



r{l)dl-— |r(Ord/ 
Jo 

(13) 

with r{l) given by Eq. (fT2|l and k{s) given by Eq. ([9]). 
Equations (fTTj) and p3|) assume that the initial condi- 
tion 2/0 is a constant. If yo is a random variable, then 
Gy [k] must be averaged over the distribution of yo to ob- 
tain the proper characteristic function of the complete 
process y{t) with its random initial condition. Here we 
shall assume that y^ — 0. 

The third and final step leading to P{Wr) consists in 
evaluating the integrals involved in Eqs. ((T2)) and p3|) . 
The integral in Eq. ([12]) defining r{l) has, with Eq. jl]), 
the solution 



r{l) 



gt;a(e('-^)/^'- - 1), 
0, 



< / < r 

1 > T. 



(14) 



From this result, we compute the first integral on the 
right-hand side of Eq. (fTS]) : 



(15) 



r(l) dl = vaq[Tr{l - e"^/^--) - r] 



As for the second integral appearing on the right-hand 
side of Eq. which involves /i, it cannot be solved an- 
alytically for all ^ S (0, 2], to the best of our knowledge, 
although it is clear that it leads to a term proportional 
to l^l^ in the exponential of Gwr {q)i since r{l) is propor- 
tional to q according to Eq. p^ . As a result, we obtain 



where 



and 



GH..(g) = c'^^«-^l«l^ 



M = v'^a[T - Tr{l - e-^/^^)], 



'0 



(16) 



(17) 



(18) 



The characteristic function of Eq. is a central 

result of this paper, which shows, following the theory 
of Levy distributions [H, [2^, , that the distribution 
P{Wr) of the work Wr is a symmetric Levy distribution 
having the following properties: 

(i) P{Wr) is symmetric and centered around M, so 
that M represents the most probable value or mode of 
Wr] see Fig.lH 

(ii) For 1 < yu < 2, the integral 



(Wr) 



W P(Wr — w) dw 



(19) 
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FIG. 1: Sketch of the distribution P{Wt) of the total (exten- 
sive) work. The mode M of the distribution is proportional 
to T while the width, which is proportional to V^^^, scales 
like r^/". 

defining the mean of Wr exists, so that M also represents 
the mean of Wr- This applies only for this range of fi 
values; for the complementary range < /i < 1, the 
integral above does not converge, so that M cannot be 
interpreted as the mean. Thus M represents the mode for 
all /i G (0, 2], but represents the mean only for /i € (1,2]. 

(iii) The parameter V is related to the width of the 
distribution P{Wr): the larger V is, the wider P{Wt) is. 
For /Lt = 2, in particular, V is half the variance of Wr, so 
the width of P{Wr), taken as the standard deviation of 
Wr, is proportional to F^^^. For < < 2, the variance 
does not exist even though V is finite. In this case, one 
can still relate V to the width of P{Wr) by calculating 
moments of Wt of order smaller than 2. In particular, it 
can be proved [30| that 

{\Wr - Mf) = aV^^", (20) 

where < (3 < fj, and a is some positive constant. The 
average (| W V — M\^) is called the fractional moment of 
order {3 [301 • ^ particular case of Eq. ([^0]) . we have 
for 13=1: 

{\Wr- M\) aV^'^ (21) 

when /i > 1. 

(iv) M and V are directly proportional to r in the limit 
of large r, that is, M ^ t and 1/ ~ r as r oo. This 
follows from Eqs.dJTI) and (HH). 

(v) For < /i < 2, P(^r) has power-law tails that 
decay according to 

(22) 

as \w\ oo [3l|. This property is responsible for the 
divergent mean observed for £ (0, 1] and the divergent 
variance observed for /x £ (0,2). Whenever, P{Wr) has 
this property, we say that the tails of P(Wr) decay as a 
power-law with exponent fj,. 

We will be interested in the next sections to study the 
fluctuations of Wt in the long-time or asymptotic limit 



where r is much larger than the relaxation time r^. To 
properly define this limit, we must note two properties 
of Wr- First, because M is proportional to r, according 
to the property (iv) above, the mode of P(Wr) escapes 
to infinity as t — > oo. Second, because V is also propor- 
tional to r, the width of P{Wt) also grows indefinitely 
with r. Thus P{Wr) gets flatter as r oo while its 
mode moves to infinity. 

To eliminate these diverging properties of P(Wt), we 
consider the intensive or scaled work, defined by Wr = 
Wr/M, as the random variable of interest rather than 
total work Wr which is extensive with r.^ The char- 
acteristic function of Wr is related to the characteristic 
function of Wr by a simple rescaling: 

GvF. (k) = Gw^ (k/M) - e'"^-^l^-l" . (23) 

With this change of variables, it is easily seen that the 
mode of Wr is M = 1, while its "width" is given by 
V = V/M'^. Therefore, the probability distribution of 
Wr is now centered at 1, and is such that V ^ r-^-f^ as 
T — > oo, since V ^ t and M ~ r in that limit. 

Note that we could have considered Wr/r as the scaled 
work rather than Wr/M. The only difference between 
these two definitions is the value of the mode M: for the 
scaled work defined by Wr/M, M is always equal to 1, 
whereas for Wr/r, M is a positive constant which may 
be different than 1. Note also that rescaling the extensive 
work Wr by r or any factor proportional to r is the only 
way of centering the distribution of Wr to a constant. 
For, if one divides Wr by a factor smaller than t, then 
the mode of the resulting random variable will still grow 
with T, whereas if one divides Wr by a factor greater 
than T, then the distribution of the resulting variable will 
become symmetric. In the latter case, all the information 
about asymmetry of the work fluctuations is discarded. 

With this in mind, we now turn to studying the dis- 
tribution of Wr for different values of /z in the range 
[0,2]. Four cases of fluctuations arising from four differ- 
ent regimes of Levy noise will be considered. The first 
is the Gaussian noise regime that leads to a conventional 
fluctuation relation for Wr- The three others are proper 
Levy noise regimes that lead to the anomalous fluctua- 
tion relation for Wr- 



B. Gaussian noise: ^ = 2 

The distribution of the scaled work Wr = Wr/M that 
results from Gaussian white noise (i.e., /i = 2) was found 
by van Zon and Cohen Q . We quickly repeat their main 
results here, since they serve as our point of departure 
for defining the anomalous fiuctuation relation. What 



^ The scaled work Wt/M is equal to Wr/iWr) when the average 
work (Wt) exists. 
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is important to note is that the distribution P(WV) for 
/i = 2 can be put in the form 



where 



P(Wr = w) 



liw) 



-r/(iu)+o(r) 



(24) 



(25) 



The paraboUc form of the function I{w), which is called 
the rate function_ , obviously leads to a Gaussian dis- 
tribution for P(Wt)- For simplicity, we shall drop in the 
remaining the o(t) term in probability distributions so 
as to write 



(26) 



P{Wr ^W)^ C-^'^'^l 



The approximation sign means, following Eq. p4)) . 
that P{Wt = w) decays, to a first degree of approxi- 
mation, exponentially with r. This property, which is 
commonly referred to as the large deviation property or 
large deviation principle (2^ , plays an important role in 
fluctuation relations, as it directly implies that 



P{Wr 



P(Wr = -W) 



where 



f{w) = I{-w) - I{w). 



(27) 



(28) 



What is special about the Gaussian noise case is that the 
fluctuation function f(w) happens to be linear in w, so 
that 



P{Wr = w) 
P(Wr ^ -w) 



(29) 



This result is the signature of the conventional fluctua- 
tion relation, as defined in [1, [l^ and Eq. ([T]) of Sec. [J 
Thus, under Gaussian noise, the mean work Wt is said 
to satisfy a conventional fluctuation relation.'^ 



C. Cauchy noise: j-i = 1 

The Fourier transform of the general characteristic 
function that we have derived for WV can be computed 
analytically for only two values of fi. The flrst value is 
(1 = 2, which we have just considered and which leads to 
Gaussian fluctuations of WV . The second value is /i = 1 



^ To be more precise, wc should say that Wt satisfies an asymp- 
totic conventional fluctuation relation, since the approximation 
above is only valid in the limit of large times t. All the results 
obtained in this paper are derived in this limit, so the attribute 
"asymptotic" will be omitted. 



and leads to a so-called Cauchy distribution having the 
form 



P(Wr = W) 



for the total work and 



V 



TT (W - M)2 -I- 1/2 



P{Wr 



V 



tt{w- Mf + 



(30) 



(31) 



for the scaled work. These distributions were already 
discussed in our previous paper [2^. What is important 
to note about P{Wr) is that its mode Af and its width 
parameter V are both proportional to t, which implies, 
as mentioned before, that this distribution moves to the 
right and flattens as t increases. This is illustrated in 
Figs.[2i;a) andllKb). The power-law tails of P{Wr) pre- 
dicted by the expression in Eq. ((22)) are clearly seen in 
the log-log plot of Fig. [2jb) as straight lines. 

The distribution P{Wt) of the scaled work retains the 
power-law tails of P(WV), as is obvious from Eqs. pO]) 
and l[3T|) as well as from Figs. [2jc) and[2l[d), but does 
not have the translation and flattening behavior that we 
have noted for P{Wr)- In fact, one exceptional property 
of P{Wr) for /i = 1 is that it becomes time-independent 
in the limit of large t because V = V/M = 0(1) in t; 
see Figs. HKc) andlDJd). This property directly implies 
that the ratio P{Wr = w)/P{Wt ~ —iB) is also time 
independent. The precise form of this ratio, which we 
denote by griw) from now on, is obtained from Eq. (j3ip : 



griw) 



P{Wr = W) {-W - 1)2 -t- 1 

P{Wr = -W) ^ (W- l)2-f 1 ' 



(32) 



and is plotted on a log-linear scale in Fig. ^e) . In this 
plot, griw) has a maximum located at iu = V2, and 



lim griw) = 1. 



(33) 



This last limit has no equivalent in the Gaussian case, 
and is a direct consequence of the fact that both tails 
of P(Wt) decay as a power- law with the same exponent 
/i = 1. In the Gaussian case, both tails of P(WV) decay 
exponentially, but because P{Wt) is not centered at 0, 
the positive fluctuations of Wt are exponentially more 
probable than the negative fluctuations. In the case of 
Cauchy noise, the difference between P{Wr = w) and 
P{Wr = —w) is so weak that the ratio griw) of these two 
probabilities goes to 1 in the limit of very large fluctua- 
tions. This implies, concretely, that large positive fiuctu- 
ations are asymptotically just as likely to be observed as 
negative fluctuations of equal magnitude, as was already 
noted in 

We shall see next that this property of Cauchy fluctu- 
ations carries over to all other values of /x in the range 
(0,2). For this reason, we follow our previous work [23[ 
and dcflne a new class of fluctuation relations based on 
this property. Given the random variable WV and its 
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FIG. 2: Work fluctuations for Cauchy noise, fj. — 1. (a) Distribution P{Wt = w) of the total (extensive) work Wt- (b) Log-log 
plot of P{Wt = w) for positive w. (c) Distribution P{Wt = w) of the scaled work Wt- (d) Log-log plot of P{Wt = w) 
for positive w. The power-law tails of P{Wt) and P{Wt) give rise to straight lines with slope —2 in the log-log plots, (e) 
Fluctuation relation: Log-linear plot of Qriw) = P{Wt ~ 'w)/P{Wt = —w). Note that P{Wt ~ w), and consequently QT-iui), 
become invariant with r for large r. This is the main feature of the Cauchy case. Units: Tt = v = b = 1. 



distribution P(WV). wc say that WV satisfies an anoma- 
lous fluctuation relation if i) P(WV) has power-law tails; 
and ii) the limit ([55)1 is satisfied. The term "anoma- 
lous" follows the terminology used in studies of Levy- type 
noise, which refer, for example, to "anomalous diffusion" 
or "anomalous transport" as diffusion or transport pro- 
cesses driven by Levy noise or noises akin to Levy noise 
(see, e.g., [12]). In this context, the term "anomalous" 
is used not in the sense of "abnormal" , but in the sense 
of "anomalous" with respect to the "normal" behavior 
obtained with Gaussian noise. 

Note that there is a further property of Gaussian fluc- 
tuations that we lose with Cauchy noise, namely, that the 
fluctuations of Wt do not decrease as r ^ oo. This again 
is a consequence of the time-independence of P{Wr) ob- 
tained with Cauchy noise. The practical consequence of 
this difference is that, in the case of Cauchy noise, one 
does not need to accumulate a large number of samples 
of Wt to observe deviations in the value of this quantity. 
For Gaussian noise, an exponential number (in r) of such 
samples is needed to observe deviations of WV from its 
mean. But for Cauchy noise any small samples of Wt 
will reveal that this quantity is fluctuating no matter 
how large t is. 



D. Upper Levy regime: 1 < /i < 2 

Two different regimes of fluctuations arise when con- 
sidering values of /i different from 1 and 2: a regime 
of "weakcr-than-Cauchy" fluctuations corresponding to 
/i e (1,2), and a regime of "stronger-than-Cauchy" fluc- 
tuations corresponding to fi S (0, 1). Wc discuss in this 
subsection the former regime, which we call the upper 
Levy regime. 

The plots shown in Fig. [3] illustrate the behavior of 
P[Wt) and P{Wt) for the case ^ = 1.5, which is repre- 
sentative of all the values fi G (1,2) in the upper Levy 
regime. These plots were obtained by numerically calcu- 
lating the inverse Fourier transform of the characteristic 
function Gwt{q): found in Eq. ((T6)) . using the "Stable" 
Mathematica package [s^. As in the Cauchy case, we see 
here that both P(Wt) and P{Wt) have power-law tails, 
and that the mode of P{Wt) increases with r, whereas 
the mode of P{Wt) is fixed at 1. However, contrary to 
the Cauchy case, the distribution P(WV) of the scaled 
work for /i € (1, 2) is not invariant with r, but becomes 
more and more concentrated around the value M = 1 
as r — > oo; see Figs. [3l[c) and[3][d). This arises because 
typical fluctuations of Wt increase less rapidly than r in 
this range of /i, which implies that the typical fluctua- 
tions for Wt must decrease with t. This concentration 
of P{Wt) around its mode implies a Law of Large Num- 
bers for Wt, in the sense that Wt 1 with probability 
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FIG. 3: Work fluctuations for jj, — 1.5 (upper Levy regime), (a) Distribution P{Wt ~ w) of the total (extensive) work Wt- (b) 
Log-log plot of P{Wt = w) for positive w. (c) Distribution P{Wt ~ w) of the scaled work Wt- (d) Log-log plot of P{Wt = w) 
for positive w. The power-law tails of P{Wt) and P{Wt) give rise to straight lines with slope —fj, — 1 in the log-log plots, (e) 
Fluctuation relation: Log-linear plot of Qriw) = P{Wt ~ w)/P{Wt ~ —w). Note that griw) — > 1 as — > oo for all t. Units: 
Tr ~ V — b = 1. 



1 as T ^ oo. This Law of Large Numbers also holds for 
Gaussian noise, and means concretely that measurements 
of Wt are most likely to be close to 1 as one considers 
longer and longer integration times r.^ 

For large but finite r, fluctuations can still be observed, 
and a measure of how positive fluctuations are more likely 
to be observed than negative fluctuations is provided by 
the ratio griw). This ratio is plotted in Fig. [3Je). As for 
the Cauchy case, we see here that positive fluctuations of 
the mean work are more likely to be observed than nega- 
tive fluctuations of equal magnitude, since Qriw) > 1 for 
w > 0, and that griw) — > 1 as — > cx), meaning that the 
difference in probabilities for positive and negative fluctu- 
ations becomes negligible for large fluctuations. Because 
of this, and the fact that P{Wr) has power-law tails, we 
conclude that Wt satisfies an anomalous fiuctuation re- 
lation for all fi G (1)2), as in the Cauchy case. Unlike 
this case, however, the maximum oi griw) increases with 
T when fi £ (lj2), and moves towards the value ui = 1 
as a result of the Law of Large Numbers. This behav- 
ior of the maximum of gr{w) is specific to /x € (1,2) 
and thus serves as a signature of the upper Levy fluc- 
tuation regime. In the Cauchy case, by comparison, the 



A random variable for which the Law of Large Numbers holds is 
also said to be "self-averaging" . 



fluctuations are stronger and lead to a time-independent 
griw) whereas, in the Gaussian case, the fluctuations are 
considerably weaker and lead to a ratio griw) which is 
exponential in r. 



E. Lower Levy regime: < ^ < 1 

The results of the numerical inverse Fourier transform 
of Gw^iq) and Gy^/ (k) are shown in Fig. 2] for /i = 0.5, 
which is representative of all the values in the range (0, 1). 
The distributions P{Wr) and P{Wr) that we obtain in 
this range of jj, characterize what we call the lower Levy 
regime, and share many of the properties that we have 
mentioned for the upper Levy regime. In particular, Wt 
satisfies an anomalous fiuctuation relation in the lower 
Levy regime, since its distribution has power-law tails 
and griw) 1 as w oo; see Fig.|Ue). 

The main difference between the lower and upper Levy 
regimes are the scaling properties of the fiuctuations with 
T. Whereas the typical fluctuations of WV decrease with 
increasing r when £ (1,2), they increase with r when 
/i e (0, 1) because, for that interval, the typical fluctu- 
ations of the "extensive" work Wr increase faster than 
r. As a result, PiWr) does not become more and more 
concentrated around its mode in the limit r ^ oo, as 
was the case in the upper Levy regime (see Fig. but 
flattens in this limit. This implies that there is no Law 
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FIG. 4: Work fluctuations for fi = 0.5 (lower Levy regime), (a) Distribution P{Wt = w) of the total (extensive) work Wt- (b) 
Log-log plot of P{Wt = w) for positive w. (c) Distribution P{Wt ~ w) of the scaled work Wt- (d) Log-log plot of P{Wt = w) 
for positive w. The power-law tails of P{Wt) and P{Wt) give rise to straight lines with slope —fi — 1 in the log-log plots, (e) 
Fluctuation relation: Log-linear plot of Qriw) = P{Wt ~ w)/P{Wr = —w). Units: Tr = v = b — 1. 



of Large Numbers for WV in the lower Levy regime: as 
longer and longer integration times are considered, the 
typical fluctuations of the scaled work actually increase 
in size, which translates in Fig. [4] into a flattening of 
griw) with T. 

Since the typical fluctuations of Wt increase with r, 
one might be tempted to rescale the total work Wr by 
a r-factor larger than r. However, the distribution that 
one would obtain from this rescaling would be symmet- 
ric: it would assign the same probability to positive and 
negative fluctuations of equal magnitude, and would thus 
contain no information about the nonequilibrium steady 
state behavior of the model, since Qt would then trivially 
be equal to 1. 



position y(t) in the comoving frame of the potential: 



IV. DISSIPATED WORK 

In previous studies of the model that we consider here, 
two quantities have been analyzed [l9j : the flrst is the 
work done by the harmonic force on the particle, which 
we have just studied; the second is the heat released by 
the particle to its environment to maintain its nonequi- 
librium steady state [11] ■ Since the particle's dynamics 
is studied in the overdamped limit, the particle has no 
kinetic energy, so that its total internal energy is entirely 
given by the potential energy, which is determined by the 



U{t) 



ix{t) - vtf y(tf 



(34) 



In this context, the heat Qt was defined in [l9| as the 
energy gained by the particle from the mechanical work 
Wt done on it by the harmonic force over a time r mi- 
nus that part of this energy which is transformed into 
potential energy, that is. 



Qt = Wt-AUt, 



(35) 



where AUt = U{t) — U{Q) is the change in potential 
energy after a time r. Here we assume that ^(0) = y(0) = 
0, so that AUt = U{t). 

It is important to note that, in the original context in 
which the model was studied [l^, the quantity Qt de- 
fined by Eq. pS]) was correctly interpreted as the heat 
because the particle is immersed in a fluid, which is re- 
sponsible for both the random force and the friction force 
applied to the Brownian particle. Accordingly, any en- 
ergy gained by the particle in the form of mechanical 
work that is not converted into potential energy is neces- 
sarily lost, by energy conservation, to the fluid as heat. In 
this sense, Eq. expresses the conservation of energy. 

In our treatment of the model, the random force ^(<) is 
an external force, which implies that the total work done 
is the sum of the work W^ done on the Brownian particle 
by the random force ^{t) and the mechanical work Wt 
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done on the particle by the moving harmonic potential. 
The heat Qr in this case should therefore be defined, by 
energy conservation, as 



(36) 



This quantity Qt is always positive because it is the heat 
produced by the friction alone, and not the friction and 
the random force as in the case of a Brownian particle 
immersed in a fluid. As a result, we cannot define a fluc- 
tuation relation for this quantity. In this section, we will 
study therefore a different quantity having both nega- 
tive and positive fluctuations. We define this quantity in 
analogy with Eq. ([35]) by 



Rr = Wr- AUr, 



(37) 



where, as before, WV is the mechanical work. The quan- 
tity Rt has a clear physical interpretation: it is that part 
of the mechanical work WV that is not converted into 
potential energy AUr- For this reason, we call Rt the 
dissipated work. Note again that Rr would be the heat 
if the only contribution to the total work done on the 
particle was the mechanical work WV; sec Eq. (f55|) . In 
our case, there are two distinct contributions to the total 
work, as shown in Eq. ([55)1 . 

The probability distribution of Wr — AUr was calcu- 
lated in the asymptotic limit by van Zon and Cohen [Toj 
for the Gaussian noise case, /i = 2. The generalization of 
their results to Levy noise is not straightforward, since 
their calculations strongly rely on the Gaussian nature 
of the noise. By assuming an independence property be- 
tween Wt and AUt for t oo, we are able, however, 
to obtain the tail behavior of P{Rr), which is sufficient 
to determine whether the dissipated work satisfies an 
anomalous fluctuation relation or not. The assumption, 
precisely, is that Wr and AUt become asymptotically 
uncorrelatcd in the limit t — > cx), so that 



-ikAUr 



) 



(38) 



in this limit. For external Gaussian noise, the distribu- 
tion of Rt obtained under this assumption can be shown, 
from the calculations reported by Taniguchi and Cohen 
[s^l, to have the same asymptotic form as the exact dis- 
tribution calculated in [1^, [l^ . We will assume here that 
this independence property between Wt and AUt holds 
also for Levy white noise when r becomes much larger 
than the relaxation time t»..^ 



^ This assumption can be argued using two basic observations. 
First, because the noise is white (delta-correlated), the position 
y(t) exhibits only short time-correlations. Second, because the 
work Wt is an integral of y{t), it can only show a weak correlation 
with any single position y{t) in time, and, in particular, with the 
last position ^(t), which determines AUt- These two arguments 
apply both to Gaussian and Levy noises. 



The characteristic function of Wt was calculated in 
Sec, mil for all /x G (0,2]. As for the characteristic func- 
tion of the potential energy. 



-ikAU-r 



) 



-^'^y'/^p{y{r)^y)dy, (39) 



it has no known closed-form solution [35}. However, it is 
known that the distribution of ?/(t) is, for large r, a Levj 
distribution with the same index /z as the noise ^{t) [3C 
so that 



Piyi^) ^v) ^ c\y\'' 



(40) 



as \y\ oo, with c a positive constant. Inserting this 
asymptotic result in the integral of Eq. (|39p . we then 
obtain 



AkAUT 



l-alfcl''/^ (41) 



as fc ^ 0. In these expressions, a and e are positive 
constants. The exact value of e is irrelevant for the last 
asymptotic result to hold, since only the tail behavior of 
P(y(r)) affects the scaling of the integral. Moreover, the 
reason for having the limit fc — > in the asymptotics of 
the characteristic function is because the behavior of the 
tails of power-law distributions is determined only by the 
behavior of their characteristic functions around fc = 0, 
and vice versa [2^. Combining this asymptotic result 
with the expression of the characteristic function of Wt 
given in Eq. p6p . we then find 

(e"=«-) ^ (1 + iMfc - V\k\^'){\ ~ a|fc|^/2) (42) 

as fc — > 0. Since /i G (0, 2), the dominant contribution in 
the above expression is 



(e"=^-) - l-a|fc|''/2^ 
which implies that 



P{Rr = r) - c' \r\ 



-1-11/2 



OO, 



(43) 



(44) 



where c' is some positive constant. Thus the tails of 
P{Rt) have a power-law decay, as in the case of the work, 
which implies that the fluctuations of the dissipated work 
satisfy an anomalous fluctuation relation similar to the 
one found for Wt- The only difference with Wt is that 
the decay exponent of the tails of P{Rt) is instead 
of ^l] sec Eq. (gH). 

This result should be contrasted with the Gaussian 
case, for which the fluctuations of Rt, or the heat Qr in 
that case, are known to satisfy an extended fluctuation 
relation p^ . [l^ , rather than a conventional fluctuation 
relation, satisfied by the work Wt- For Levy noise, the 
fluctuations of Wt are not so different from the fluctua- 
tions of Rt , since both quantities have distributions with 
power-law tails. In spite of this difference, there is one 
property of Rt that remains the same for both Gaussian 
and Levy noises, namely, that the large fluctuations of 
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Rt are mostly the result of the large fluctuations of the 
potential energvALL. This property, also observed for 
Gaussian noise [la, [3, can be understood here by noting 
that the asymptotic behavior of the characteristic func- 
tion of Rr around the origin fc = 0, which determines 
the behavior of the tails of P{Rr), is determined entirely 
by the asymptotics of the characteristic function of AUr 
around fc = 0. Thus, although Rr has different fluctua- 
tion properties depending on the noise used (Gaussian or 
Levy) , its fluctuations are mostly the result of the fluc- 
tuations of AUt no matter what noise is applied. 

V. FLUCTUATION RELATIONS AND LARGE 
DEVIATIONS 

It should be clear from the previous results that what 
differentiates conventional and extended fluctuation re- 
lations, on the one hand, and anomalous fluctuation re- 
lations, on the other hand, is the existence of a large 
deviation principle for the distribution of the observable 
studied. In the case of the scaled work Wr, for example, 
P{Wr) satisfies a large deviation principle for Gaussian 
noise, as reported in Eq. (j24|) . but not for strict Levy 
noise with /i G (0,2). In the case of Cauchy noise, in 
particular, the distribution P{Wt) is such that the limit 

/(w) = lim --\nP{Wr ^w) (45) 

yields I(w) = for all w, which means in effect that 
P{Wr) does not satisfy a large deviation principle. The 
same result applies to the distribution of WV obtained 
for all fi G (0,2), as well as for the distribution of the 
dissipated work R^. obtained for all fj, € (0,2), because 
all these distributions have power-law tails in this range 
of Id. This result is also general insofar as any observable 
Ar that does not obey a large deviation principle does not 
obey a conventional or an extended fluctuation relation, 
as is obvious from the definition of these two types of 
fluctuation relations given in Eqs. ([1]) and ([2]). Whether 
or not At satisfies in this case an anomalous fluctuation 
relation depends on the explicit form of P{Ar), but we 
know for sure that A.,- satisfies neither a conventional nor 
an extended fluctuation relation. 

This last observation brings us to the question of 
whether there exist fluctuation relations other than con- 
ventional, extended, and anomalous. In other words, if 
an observable At does not satisfy a conventional or an 
extended fluctuation relation, does it necessarily satisfy 
an anomalous fluctuation relation? 

It is difficult at this point to answer this question, since 
there is nothing in principle that prevents one from imag- 
ining noises that lead to distributions that are not ex- 
actly Levy and yet do not obey a large deviation princi- 
ple. However, the fact that there exists a link between 
fluctuation relations and large deviation theory restricts 
somehow what can be imagined. Indeed, it is known 
from this theory that random variables obeying large de- 



viation principles have, in most cases, flnite moments at 
all order. Therefore, it is natural to expect that noises 
having finite moments should lead to conventional or ex- 
tended fluctuation relations, because they should lead to 
distributions having a large deviation form. This is the 
case for Gaussian noise, as we have seen here, but also 
for other noises having finite moments, including Pois- 
son noise; see [36i] . On the other hand, noises having 
one or more infinite moments should lead to anomalous 
fluctuation relations, since they cannot lead in general to 
distributions having a large deviation principle. This is 
the case for Levy noise, and it should be the case, too, 
for any noises having distributions with power-law tails. 
From this point of view, the results that we have obtained 
for Levy noise should be representative of a larger class of 
fluctuations arising from any noises having one or more 
infinite moments. 



VI. DISCUSSION 

We close this paper with some remarks about possible 
experimental verifications of our results, some technical 
issues about Levy noise, and extensions of our results to 
nonlinear models and other power-law noises. 

(1) The fact that Levy noise has an infinite noise power 
implies, as we have mentioned in Sec. |TT1 that one can- 
not define a fluctuation-dissipation relation between the 
friction force in the model considered here and the power 
of the noise. But, as argued in that section, this is not 
a problem insofar as the noise is external, i.e., that it is 
produced and imposed externally by a physical process 
which is different from the physical process giving rise to 
the friction. Any experimental verification of our results 
will have to include a noise source that has this property 
of being external. Perhaps the easiest way to produce 
such a noise is to generate Levy white noise on a com- 
puter (see the remark (4) below) , and to feed the output 
to an appropriate transducer (electrical or mechanical), 
which will transform the digital noise produced by the 
computer into a physical noise, e.g., an electric noise or 
a mechanical noise. 

(2) We have proposed in our previous paper [2^ two 
experiments that could be used to test our results. These 
experiments follow the suggestion of the previous para- 
graph that the noise must be external. 

In the first experiment, a solid object is placed on a 
solid table that is vibrated horizontally by a mechanical 
transducer controlled by a computer; see [33]. By vibrat- 
ing the table with Gaussian noise or by Levy noise, one 
should be able, in principle, to generate a steady-state 
motion of the object, and to study the fluctuations of 
work- and heat-related quantities, such as the work done 
on the object by the gravitational force in the case where 
the table is slightly tilted. An experiment of this sort is 
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under development.^ The requirement that the friction 
and the noise are uncoupled is obviously satisfied in this 
case, since the friction is the solid-solid friction, which 
has nothing to do with the arbitrary noise externally im- 
posed by the vibrating table. 

The second experiment that we proposed in [23| is 
based on granular gases maintained in steady states by 
vibrating their container. Many studies have looked at 
the properties of these gases when they are vibrated 
by periodic for cing and by Gaussian noise (see, e.g., 
[ill . [T^ . [Tsl . [l^ . [ssfl. A natural variation of these ex- 
periments, in view of our work, is to change the forcing 
signal by Levy white noise. This can be done, in prin- 
ciple, in experiments, and certainly in numerical simu- 
lations of granular gases. The quantity that is usually 
studied for these systems is the power injected by the 
vibrating force. As for the mechanical work studied in 
this paper, one could study the power injected, and ver- 
ify that the fluctuations of the latter quantity satisfy an 
anomalous fluctuation relation for Levy noise. 

(3) A different experiment, which could be used to 
probe the Levy fluctuations of i?r, can be imagined us- 
ing the analogy that exists between the fluctuations of 
the Brownian particle studied in this paper and the cur- 
rent fluctuations of small RC electrical circuits @, [lol |. 
The noise in such circuits is internal, since it is the ther- 
mal Johnson-Nyquist noise, usually treated as Gaussian 
white noise. However, there is nothing that prevents 
one from introducing an additional, external noise source 
in these circuits via fluctuating voltage sources, such as 
those studied, e.g., in [s^. The fluctuations of the volt- 
age sources can be generated by a computer, which means 
that they can be used, in principle, to mimic Levy white 
noise. Under Gaussian (thermal Johnson-Nyquist) noise, 
it is known that the fluctuations of the heat dissipated by 
the resistance in RC circuits follow an extended fluctua- 
tion relation [l,[l3|- Based on our results, we expect that, 
in the presence of an additional Levy noise, the heat fluc- 
tuations will be similar to those of the dissipated work 
Rr studied here. 

(4) For the purpose of the experiments just described. 
Levy noise can be generated physically to a good degree 
of accuracy by generating a Levy noise on a computer, us- 
ing techniques similar to those used for generating Gaus- 
sian noise (see, e.g., [23, HO; El), and by feeding this 
artificial noise into a mechanical or electrical transducer. 
Of course, Levy noise generated on a computer is never 
exactly "Levy": like any physical noise, there is always 
a maximum value of the noise that can be applied to a 
system. In the case of Levy noise, this maximum value or 
cut-off transforms the noise into a truncated Levy noise 
[4^ for which the distributions of WV and Rr do not have 
the exact form of a Levy distribution, simply because 
truncated Levy noise has finite moments to all order. 




M. K. Chaudhury, private communication, 2009. 



FIG. 5: Convergence of distributions for truncated Levy noise. 
(Colored lines) Log-log plot of P{Wt) for truncated Cauchy 
noise with various cut-off values ^max, which is the maxi- 
mum value imposed to the noise ^{t). The curves are ob- 
tained by a direct numerical sampling of the trajectories of the 
Langevin equation ^ with truncated Cauchy noise. (Black 
line) P{Wt) obtained with exact Cauchy noise, as calculated 
in Eq. (|3ip . Units: Tr = v = b — 1. Simulation parameters: 
At — 10^^ sec, r = 5 sec. 



However, it is possible to reveal the Levy character of 
these distributions by studying their behavior or "trend" 
as the cut-off is increased. 

Figure [S] shows, as an illustration, the distribution of 
Wt obtained by direct sampling of Eq. ^ using trun- 
cated Cauchy noise with different cut-off values. The 
behavior of P(Wt) seen in this figure is general for trun- 
cated Levy noise [i^jH^I' i-G-, as the cut-off is increased, 
the distribution of P(WV) approaches the Levy distribu- 
tion predicted for exact Levy noise, which in this case 
is the Cauchy distribution given in Eq. ([3T|) . The same 
behavior is expected to arise for Rr- 

(5) For strict Levy noise with fi E (0,2), there is al- 
ways a very large fluctuation that dominates the other 
fluctuations in time, especially in the lower Levy regime, 
/i G (0, 1). For Gaussian noise, on the other hand, all ob- 
servable fluctuations are more or less of the same order of 
magnitude. This difference between Levy and Gaussian 
noise should directly be observable in the experiments 
mentioned above. 

(6) The generalization of our results to nonlinear 
Langevin equations with Levy noise should give rise to 
interesting results. It is known, for example, that the 
stationary distribution of a Langevin equation involving, 
as here, a quadratic potential is unimodal for Gaussian 
and Levy noise. In the case of a quartic potential, how- 
ever, the stationary distribution of the position is still 
unimodal for Gaussian noise but bimodal for Levy noise 
[4^ . An interesting question, in the context of this re- 
sult, is whether the distributions of WV and Rr obtained 
for a moving quartic potential are also bimodal. For the 
quartic potential, it is also known that the distribution 
of the position has a finite variance, so it would be inter- 
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esting to see whether or not the fluctuations of WV and 
Rt are anomalous in this case. 

(7) A recent paper by Chechkin and Klagcs has ap- 
peared recently, which studies an anomalous fluctuation 
relation similar to the one studied in this paper and our 
previous work [l^. The paper by Chechkin and Klages 
deals with Levy noise, as in this paper, but also with 
long-time correlated Gaussian noise, which gives rise to 
a probability distribution of the work having a stretched- 
exponential form [isj . 
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